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(<— ^ I Using the asymptotic-expansion technique, we compute the dominant two-loop 

electroweak corrections, of O^Gprrif), to production and decay via a pair of pho- 
tons or gluons of the CP-odd Higgs boson A^ in a two-Higgs-doublet model with 
r> I low- to intermediate values of the Higgs-boson masses and ratio tan/3 = ^2/^1 of 

C^ ' the vacuum expectation values. We also study the influence of a sequential heavy- 

fermion generation. The appearance of three 75 matrices in closed fermion loops 
requires special care in the dimensional regularisation of ultraviolet divergences. 
The finite renormalisation constant for the pseudoscalar current effectively restor- 
ing the anticommutativity of the 75 matrix, familiar from perturbative quantum 
chromodynamics, is found not to receive a correction in this order. We also revisit 
the dominant two-loop electroweak correction to the /7 — > 77 decay width in the 
standard model with a fourth fermion generation. 
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1 Introduction 

The search for Higgs bosons is among the prime tasks at the Fermilab Tevatron and will 
be so at the CERN Large Hadron Collider (LHC), to go into operation later during this 
year, and the International e"'"e~ Linear Collider (ILC), which is currently being designed. 
The standard model (SM) contains one complex Higgs doublet, from which one neutral 
CP-even Higgs boson [H) emerges in the physical particle spectrum after the electroweak 
symmetry breaking. Despite its enormous success in describing almost all experimental 
particle physics data available today, the SM is widely believed to be an effective field 
theory, valid only at presently accessible energy scales, mainly because of the naturalness 
problem related to the fine-tuning of the cut-off scale appearing quadratically in the Higgs- 
boson mass counterterm, the failure of gauge coupling unification, the absence of a concept 
to incorporate gravity, and the lack of a cold-dark- matter candidate. Supersymmetry 
(SUSY), which postulates the existence of a partner, with spin shifted by half a unit, to 
each of the established matter and exchange particles, is commonly viewed as the most 
attractive extension of the SM solving all these problems. The Higgs sector of the minimal 
SUSY extension of the SM (MSSM) consists of a two-Higgs-doublet model (2HDM) and 
accommodates five physical Higgs bosons: the neutral CP-even h^ and H^ bosons, the 
neutral CP-odd A^ boson, and the charged ilf^-boson pair. At the tree level, the MSSM 
Higgs sector has two free parameters, which are usually taken to be the mass M^o of the 
A^ boson and the ratio tan/5 = ^2/^1 of the vacuum expectation values of the two Higgs 
doublets. 

The discovery of the A^ boson would rule out the SM and, at the same time, give 
strong support to the MSSM. At the LHC, this will be feasible except in the wedge of 
parameter space with M40 ^ 250 GeV and moderate value of tan /?, where only the /i° 
boson can be detected p]. For low to intermediate values of tan/?, gluon fusion is by 
far the dominant hadroproduction mechanism. At large values of tan/5, A'^bb associated 
production becomes important, too, especially at LHC cm. energy, y/s = 14 TeV ^^3\■ At 
the ILC operated in the 77 mode, via Compton back-scattering of highly energetic laser 
light off the electron and positron beams, single production of the A^ boson will allow 
for its discovery, also throughout a large fraction of the LHC wedge, and for a precision 
determination of its profile [4] . Two-photon collisions, albeit with less luminosity, will also 
take place in the regular e^e~ mode of the ILC through electromagnetic bremsstrahlung 
or beamstrahlung off the lepton beams. 

In the mass range M^o < 2mt and for large values of tan/3 in the whole M^o range, 
the A^ boson dominantly decays to a bb pair, with a branching fraction of about 90% 
[2][3] . As in the case of the H boson of the SM, the rare 77 decay channel then provides 
a useful signature at the LHC if the b and b quarks cannot be separated sufficiently well 
from the overwhelming background from quantum chromo dynamics (QCD). The A^ -^ gg 
channel will greatly contribute to the decay mode to a light-hadron dijet, which will be 
measurable at the ILC. 

Since the A^ boson is neutral and colourless, the ^4^77 and A^gg couplings are loop 
induced. As the A^ boson has no tree-level coupling to the W boson and its coupling to 



sfermions flips their "handedness" (left or right), the ^^77 coupling is mediated at leading 
order (LO) by heavy quarks and charged leptons and by light charginos, whereas heavy 
charginos decouple |5j. The A'^gg coupling is generated at LO by heavy-quark loops [6]. 

Reliable theoretical predictions for the y4''77 and A^gg couplings, including higher- 
order radiative corrections, are urgently required to match the high precision to be reached 
by the LHC and ILC experiments [71I8] . Specifically, the properties of the A^ boson, 
especially its CP-odd nature, must be established, and the sensitivity to novel high-mass 
particles circulating in the loops must be optimised. The present state of the art is 
as follows. The next-to-leading-order (NLO) QCD corrections, of relative order 0{as) 
in the strong-couphng constant ««, to the partial decay widths T{A'^ -^ 77) [9p0] and 
r(A° — > gg) [To], and the production cross section a{gg — * A^) [TOfTT] are available for 
arbitrary values of quark and A°-boson masses as one- dimensional integrals, which were 
solved in terms of harmonic polylogarithms for r{A^ -^ 77), r{A^ -^ gg), and the virtual 
correction to cr{gg -^ A^) |12i,13j . The latter was also obtained for general colour factors 
of the gauge group SU(iVc) in the limit m^ — * 00 using an effective Lagrangian [T3|. In 
the same way, the 0{as) correction to a{gg -^ A^) was first calculated in Ref. [T5] . 

The next-to-next-to-leading-order (NNLO) QCD corrections, of (9(a^), to r{A^ —>■ 
gg) [IS] and (y{gg — > A^) [T^ were found for m^ — >■ 00 using an effective Lagrangian. 
The 0{as) SUSY QCD correction, due to virtual squarks and gluinos besides the heavy 
quarks, to a{gg -^ A°) was obtained from an effective Lagrangian constructed by also 
integrating out the SUSY particles [18]. The two- loop master integrals appearing in the 
latter calculation if the masses of the virtual scalar bosons and fermions are kept finite 
were expressed in terms of harmonic polylogarithms ^^ . 

In this paper, we take the next step and present the dominant electroweak corrections 
to r(A° -^ 77) and r(A° — >• gg) at NLO. Our key results were already summarised 
in Ref. [T^. Here, we present the full details of our calculation and a comprehensive 
discussion of its phenomenological implications. Since these corrections are purely virtual, 
arising from two-loop diagrams, they carry over to (7(77 —>■ A^) and (y{gg -^ A^), via 

a(77M - ^°) = j;^?^nA' - ii/gg)6 {s - Mlo) , (1) 

where A^^ = 1 and A^^ = A^^ — 1 = 8 are the colour multiplicities of the photon and 
the gluon, respectively, and s is the partonic cm. energy square. For the time being, 
we focus our attention on the particularly interesting region of parameter space with low 
to intermediate Higgs-boson masses, M/jO, M/^o, M^o, M//± < mtjj and low to moderate 
value of tan/5, tan/? -C mt/mb, and assume that the SUSY particles are so heavy that 
they can be regarded as decoupled, yielding subdominant contributions. The dominant 
electroweak two-loop corrections are then induced by the top quark and are of relative 
order 0{xt), where Xt = G pm'f / (Stt"^ \/2) ^ 3.17 x 10"^ with Gp being Fermi's constant. 
We also consider the influence of a sequential generation of heavy fermions F, beyond the 



^As for Mj^o, we actually need M^o < 2Afiy±,2Af^± in order for asymptotic expansion to be appli- 
cable. 



established three generations, which generate corrections of generic order 0{xf)- Such 
corrections were aheady studied for the if — i> 77 decay in the SM supplemented with a 
fourth fermion generation in Ref. [20], and we revisit this analysis. 

In the calculation of two-loop electroweak corrections to the 74*^77 and A^gg couplings, 
one encounters closed fermion loops involving three 75 matrices, so that the use of the naive 
anticommuting definition of the 75 matrix is bound to fail. This leads us to employ the 
't Hooft-Veltman-Breitenlohner-Maison (HVBM) [21] scheme and a finite renormalisation 
constant, Z^, for the pseudoscalar current to effectively restore the anticommutativity of 
the 75 matrix [22|23|l24j . which is so far only known within QCD [2^ . 

This paper is organised as follows. In Section [21 we explain our method of calculation 
and evaluate Z^ to 0{xt) and 0{xf)- In Section [3l we calculate the 0{xt) and 0{xf) 
corrections to r(A° -^ 77) and r(y4° -^ gg) in the 2HDM with three and four fermion 
generations. We also recover the 0{as) correction to r(A° — >■ 77). In Section HI we 
recalculate the 0{xf) correction to r(if — >• 77) due to a fourth fermion generation added 
on top of the SM. In Section we present our numerical results. We conclude with a 
summary in Section [6l 

2 Method of calculation 

As explained below, we assume a strong hierarchy between the heavy fermions on the one 
hand and the gauge and Higgs bosons on the other hand, so that we may extract the 
leading corrections using the asymptotic-expansion technique [25]. We use a completely 
automated set-up, which relies on the successive use of the computer programs QGRAF [26] , 
q2e, exp [27], and MATAD [2H]- First, QGRAF is used to generate the Feynman diagrams. 
Its output is then rewritten by q2e to be understandable by exp. The latter performs the 
asymptotic expansion and generates the relevant subdiagrams according to the rules of 
the so-called hard- mass procedure [25]. Form [29] files are generated. They can be read 
by MATAD, which performs the calculation of the diagrams. 

Since we consider the SUSY partners to be decoupled, we may as well work in the 
2HDM without SUSY. We may thus extract the ultraviolet (UV) divergences by means 
of dimensional regularisation, with d = 4 — 2e space-time dimensions and 't Hooft mass 
scale n, without introducing SUSY-restoring counterterms [30]. For convenience, we work 
in 't Hooft-Feynman gauge. We take the Cabibbo-Kobayashi-Maskawa quark mixing 
matrix to be unity, which is well justified because the third quark generation is, to good 
approximation, decoupled from the first two [31]. We adopt Sirlin's formulation of the 
electroweak on-shell renormalisation scheme [53|, which uses Gp and the physical particle 
masses as basic parameters, and its extension to the MSSM [33]. Various prescriptions 
for the renormalisation of the auxiliary variable tan/3, with specific virtues and flaws, 
may be found in the literature, none of which is satisfactory in all respects (for a review, 
see Ref. [S])- For definiteness, we employ the Dabelstein-Chankowski-Pokorski-Rosiek 
(DCPR) scheme [2SIM], which maintains the relation tan/3 = ^2/^1 in terms of the 
"true" vacua through the condition 6vi/vi = SV2/V2, and demands the residue condition 



ReE^o(Af^o) = and the vanishing of the A^-Z'^ mixing on shell as ReS^()2()(M^o) = 0, 
where S^o(g^) and S^02o(g^) are the renormalised y4*^-boson self-energy and A'^-Z^ mixing 
amplitude, respectively. It has been pointed out [31] that the DCPR definition of tan/3 
is gauge dependent. We do not actually encounter this drawback in our analysis, since 
we need to renormalise tan/5 to 0{xf), so that only fermion loops contribute. However, 
this problem will show up when subleading terms of the two-loop electroweak corrections 
are to be computed. Our final results can be straightforwardly converted to any other 
renormalisation prescription for tan/5, by substituting the finite relationship between the 
old and new definitions of tan/?. 

As already mentioned in Section [H the evaluation of the relevant two-loop diagrams 
is aggravated by the appearance of three 75 matrices inside closed fermion loops. This 
leads us to adopt the HVBM scheme [21], which allows for a consistent treatment of the 
Dirac algebra within the framework of dimensional regularisation. In this scheme, the 75 
matrix is given by 

75 = ^e/..p.7V7V, (2) 

where the totally antisymmetric Levi-Civita tensor is defined in d dimensions as 

if (/i, z/, p, a) even permutations of (0,1,2,3), 
^fiupa- = ■{ ~1 if (/^)^)P)Cr) odd permutations of (0,1,2,3), (3) 

otherwise. 

This definition leads to the following mixed anticommutation and commutation relations, 
where we have to distinguish between 4 and {d — 4) dimensions: 

{7^75} = 0, if /i = 0,1, 2, 3, 
[7^5 75] = 0, otherwise. (4) 

By giving up the full anticommutation property of 75, we can retain the familiar expression 
for the trace of four 7 matrices and one 75 matrix: 

tr(7^7!.7p7a75) = 4ie^^p^. (5) 

Traces involving less than four 7 matrices and one 75 matrix vanish. We introduce the 
following projectors onto the 4- and {d — 4)-dimensional subspaces: 



r" 



g^'^ if fx and u are smaller than 4, 
otherwise; 



-;ji/ I 9^'^ if /^ ^^^ ^ ^^6 larger than 3, 

9 = \ (6) 

I otherwise. 

Here and in the following, we label quantities in 4 dimensions with a tilde, quantities in 
{d — 4) dimensions with a hat, and quantities in d dimensions without superscript. For a 



given four- vector V, we thus have 



yt^ = ^/^^v;. (7) 

The projectors fulfil the following relations: 

r, = g'^-g,. = r~g,. = 4, 
g', = g^''g,. = g^''g,u = d-4, 

r'g.u = 0. (8) 

More details may be found in Ref. [23j. We explicitly verified that, in our case, the naive 
anticommuting definition of the 75 matrix yields ambiguous results, which depend on the 
way of executing the Dirac traces. 

The actual implementation of these rules into the MATAD setup is accomplished in two 
independent ways. Firstly, we use the Mathematica package TRACER [37j and compute 
the traces separately. Secondly, we implement our own FORM routine for evaluating the 
traces in the HVBM scheme. Both methods yield the same results. 

The application of the HVBM scheme introduces loop momenta that are projected onto 
the {d — 4)-dimensional subspace. These have to be expressed through loop momenta in 
the full d-dimensional space because MATAD performs the integration in d dimensions. For 
instance, in the case of a massive one-loop tadpole with loop momentum q, we substitute 

f = -/ (9) 

in the numerator of the integrand. Similar identities can be derived for all the other cases. 
Furthermore, we have to introduce an additional finite counterterm, Z5, in the renor- 
malisation of the pseudoscalar current to effectively restore the anticommutativity of the 
75 matrix P^23f24j . Within QCD, Zf is known through 0(a^) [21]. Here, we need Z5 
at 0{xt) and 0{xf)- In order to explain our procedure, we first repeat the derivation of 
the 0{as) term. For simplicity, we work in massless QCD with 't Hooft-Feynman gauge. 
As usual, we adopt the modified minimal-subtraction (MS) renormalisation scheme. The 
pseudoscalar current is defined in coordinate space as 

P{x) = Z2ZPZlij{xh,ij{x), (10) 

where Z2 and Z^ are the usual wave-function and pseudoscalar-current renormalisation 
constants, respectively. Passing to momentum space, we define 

iSMT,{p,p')iSFip') = [d''xd%e'^P-^-P'y\0\T[^P{x)PiOmy)]\0), (11) 



where T denotes the time-ordered product and iSpip) = J d'^xe^P'^{0\T['ip(x)'i(j{0)]\0) is 
the Feynman propagator of the quark. The key quantity for our purposes is then the 

6 



amputated Green function T^lpjp) at zero momentum transfer. The Feynman diagrams 
relevant through 0{as) are depicted in Fig. [H where crosses and dots indicate the inser- 
tions of P{x) and the operator renormahsation, respectively, and it is understood that 
external legs are amputated. The tree-level diagram in Fig. [TJ^a) yields 

tP{p,p)=j,. (12) 

The one-loop diagram in Fig. [U^b) leads to the integral 

r, {p,p) = 4na.C,J (^(^3^p^, (13) 

where Cp = {N^ — l)/{2Nc) = 4/3 is the eigenvalue of the Casimir operator of the 
fundamental representation of SU(A^c), with Nc = 3 for QCD. We decompose the string 
of gamma matrices in the numerator as 

The first term on the right-hand side of Eq. flT4l) is what we would obtain using an 
anticommuting 75 matrix. Upon loop integration it yields an expression proportional to 
Eq. flT^ , the divergent part of which is 

Ti^'''\p.p) = -C,\. (15) 

This is removed by Fig. [11(c), the 0{as) terms of Z2 and Z^ being [Mj 

Zo ^1 — — G p — , 
TT 4e' 

ZP = l- —Cf — . (16) 

TT 4e 

The remaining terms on the right-hand side of Eq. (IT^ are evanescent; they live in the 
unphysical {d — 4)-dimensional part of space-time and vanish if we let d — * 4. Yet the 
loop integral is divergent, so that an unphysical finite contribution remains. We can 
ensure the vanishing of such contributions to all orders by a finite renormahsation. This 
is exactly what is achieved by the finite renormahsation Zf. In this way, we ensure that 
the evanescent part does not mix into the physical part. Upon loop integration, the 
evanescent part of Eq. flT^ yields 

r«'-^(p,p) = 2^Cp75, (17) 

TT 

SO that, through C(as), the finite counterterm is 

ZI = 1 + 5ZI=1-2'^Cf, (18) 

in agreement with Ref. [2] 



(a) (b) (c) 

Figure 1: Feynman diagrams contributing to Z^ at ©(os). Crosses and dots indicate insertions 
of P{x) and its operator renormalisation ^2^^^!*, respectively. 



We now apply the same procedure at Oixt) and Oixp). To this end, we have to 
consider the counterparts of the diagram in Fig. [I](b) where the gluon is replaced by 
neutral or charged scalar electroweak bosons, S^ = ■}^,h^,H^,A^ and S"^ = 0^,iif^, 
where x° and 0^ denote the Goldstone bosons. These are depicted in Figs. [2](a) and 
(b), respectively. We find that the sets of diagrams in Figs. [2](a) and (b) add up to zero 
separately. In the three-generation case, the diagrams in Figs. [2](b) do not contribute in 
0{xt) at all. Consequently, we have 6Z^ = at 0{xt) and 0{xf). 

/ / / \ 

(a) (b) 

Figure 2: Feynman diagrams contributing to Z^ at 0{xt) and 0{xf). Crosses indicate insertions 
of P{x), and 5° = x°, /i°, H^, ^°, S^ = </>±, H^, and f = t,b, U, D, N, E denote generic neutral 
and charged scalar bosons and fermions, respectively. 
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3 A^ ^- 77 and A^ -^ gg 



In this section, we first discuss the case of A^ -^ 77 in detail. Specifically, we recall the 
Born result and recover the 0{as) correction in Subsection 13.11 and evaluate the 0(xj) 
and 0{xf) corrections in Subsection 13.21 In Subsection 13.31 wg then extract from the 
latter the corresponding corrections for A^ -^ gg. 

We drop all terms including an even number of 75 matrices in the fermion trace, since 
they do not give contributions to the final results. Those with an odd number lead to 
expressions that are proportional to the epsilon tensor. The transition amplitude can thus 
be decomposed as follows: 

r = ya,,se;{qi)e:{q2)A^-^'^^^', (19) 



where qi and e*(gi) are the four-momenta and polarisation four- vectors of the outgoing 
photons i = 1,2. By Lorentz covariance, Eq. (fT9l) can be written as 



r 



_flUp(T 



<lif^(l2ue*p{qi)el{q2)A, 



(20) 



^P^a 



where 

/] _ 9[ati9(3ugypgScr]iC['lC[2 Afiuaf3-y5 /^i \ 

48(gi-g2)^ 
Here, [ ] i denotes antisymmetrisation in the first indices of the metric tensors, and it is 
understood that the external momenta gi and g2 have non-vanishing components only in 
the physical four dimensions of space-time. The partial width of the A^ -^ 'j'j decay is 
then given by 

r(A°-77) = ^l^r. (22) 



647r ' ' 
The form factor A is evaluated in perturbation theory as 



-4 = E(-^ 



LO 



Af 



X/ + . ..)+..., 



(23) 



where the sum is over heavy fermions f = t,F, Ah'^ is the one-loop contribution induced 
by charged fermions, A"^" is the two-loop QCD correction in the case of / being a quark, 
A/ is the dominant two-loop electroweak correction due to weak-isospin doublets of 
quarks and leptons, and the ellipses stand for the residual one- and two-loop contributions 
as well as all contributions beyond two loops. The counterpart of Eq. (!22|) for T{A^ -^ gg) 
contains an additional colour factor of Ng/A = 2 on the right-hand side. 

The couplings of the A^ boson to fermions are proportional to their masses. Therefore, 
we only consider diagrams where the A^ boson couples to the top quark or a fourth- 
generation fermion. However, we must bear in mind that diagrams where it couples to 
the bottom quark may become sizeable for large values of tan/3, for tan/5 = Olmt/mh), 
because that coupling is proportional to tan/5. In the following, we thus concentrate on 
low to intermediate values of tan/?. 



3.1 Born result for T{A^ — ^ 77) and 0{cXs) correction 

The LO result arises from the diagrams shown in Fig. [3] and may be found in Ref. [5] for 
d = 4 space-time dimensions. For d = 4 — 2e through 0{e), we have in closed form and 
as an expansion in Tf = M^o/(2m/)^: 



J ^ 



71 



)9f 



)9f 



47r/i2 



m. 



-JE 



A-njji^ 



lE 



m. 



y 



arcsin^ ^/tJ + 0(e) 
1 8 2 4 3 128 4 



1575 



(24) 




Figure 3: One loop diagrams contributing to T{A^ -^ 77). f = t,b, U,D,N,E denotes generic 
fermions. 



where je is the Euler-Mascheroni constant, a^m is Sommerfeld's fine-structure constant, 
Nf = l,Nc is the colour multiphcity of ferniion /, Qj is its fractional electric charge, 
and gj = cot/?, tan/? for up-type and down-type is its coupling strength to the A^ boson 
normalised to its Yukawa coupling in the SM. Note that M^o enters Eq. i^^ through the 
kinematic relation (gi + ^2)^ = M'^g. The corresponding result for A^ -^ gg is obtained 
from Eq. flM|) through the substitution aemNtQ'l — > a^. 




.g A^ q 




g q 



Figure 4: Sample two-loop diagrams contributing to the 0{as) QCD corrections to r(^'' — > 77). 



We now turn to the 0{as) QCD corrections. Besides the proper vertex diagrams, some 
of which are depicted in Fig. HJ we also need to include the counterterm contribution, so 
that Ag'' = ^qCT + -^9,0- ^s for renormalisation, besides including the finite counterterm 
of the pseudoscalar current in Eq. fITH]) . we only need to renormalise the quark mass 
appearing in the prefactor of Eq. fl24|) . by shifting its bare value as m° = rriq + Srriq, so 
that 

^LO| Jsz'',-2e^-^]. (25) 






m„ 



In the on-shell scheme, the quark mass counterterm is 



6mn 



m„ 



Of 

TT 



3 ^ 3 , /i^ 



mt 



(26) 
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where A = 1/e — 7s + ln(47r). In total, our evaluation yields 



vr TT 



16 



68 2 _ 53012 3 

45^" ^ 42525^" 



34712 
33075''^ 



rt + (^(^,') 



(27) 



which agrees with the Taylor expansion of the analytic result derived in Refs. [T^IITB] from 
the integral representation originally obtained in Refs. PITU] . 

Notice that the 0{t^) term in Eq. (1271) vanishes, so that the 0{as) correction is 
suppressed for small values of M^o. In fact, as a consequence of the Adler-Bardeen 
theorem [38] , the large-mt effective Lagrangian of the ^''77 interaction does not receive 
QCD corrections at any order pil5f39j . 



3.2 0{xt) and 0{xf) corrections to r(A° — ^ 77) 




A' 



f 




f 



f 



7 



S 



7 






Figure 5: Generic two-loop diagrams contributing to the 0(xf) electroweak corrections to 
T{A^ -^ 77). S = x°,<t>^,h'^,H°,A°,H^ and / = t,b,U,D,N,E denote generic scalar bosons 
and fermions, respectively. The couplings of the neutral scalar bosons to the bottom quark are 
to be neglected and those of the neutral particles to the photon vanish. 



Now we turn to the dominant electroweak two-loop corrections. We first consider 
the case of three fermion generations. Later on, we also study the additional corrections 



11 



induced by a sequential fermion generation, consisting of an up-type quark U, a down- 
type quark D, a Dirac neutrino N, and a charged lepton E. The contributing diagrams 
are depicted generically in Fig. [5l We have checked exphcitly that diagrams including 
virtual W bosons do not contribute to our order. 

Let us first discuss the counterterm contributions. In contrast to the QCD case in 
Eq. f l25|) . we now also have to renormalise the A^ wave function, Gp, and tan/3. These 
additional contributions are universal. We thus have 

IT 



4^CT = -2'^'G'J'^NfQ}gf (^6ZI - 2e^ + <5.) 



:28i 



where, in the electroweak on-shell scheme [321133] supplemented with the DCPR [35] def- 
inition of tan/3. 

Here, 5v/v is the common DCPR counterterm for the two Higgs doublets given in Eq. (3.11) 
of Ref. [36] and Ar [32] contains those radiative corrections to the muon lifetime which 
the SM introduces on top of those derived in the QED-improved Fermi model. In terms 
of (transverse) self-energies, we have 



Su 



Siy±,T(0) 



- J:'^oiMlo) + (tan/3 - cot^) 



I!^o^o(M^o) 



(30) 



In the three-generation case, we set rrih 
hierarchies: 



and formally impose the following mass 



Mz,Mw, MhO , M^o , M^o , Mh± < rrit , M^o <2Mw, 2Mh± , 



(31) 



which ensure the applicability of the asymptotic-expansion technique. In practice, this 
implies that the unknown Higgs-boson masses obey 



Mho,MHO,MAO,MH± <mt, 



Mao < 2Mh±. 



(32) 



The leading two-loop electroweak corrections are then of 0{xt). In the presence of fourth- 
generation fermions F = U,D, N, E, we assume that their masses obey 



Mz, Mw, MhO,MHO, Mao, Mh± < mu, mo, m^, m^. 



(33) 



For simplicity, we consider the special cases mu 3> rriD, mu = m/5, and niu <C mu, and 
similarly for the leptonic weak-isospin doublet, so that we are effectively dealing with 
single-scale problems yielding corrections of 0{xf)- 

We first list the non-vanishing counterterms entering Eq. ( 128|) for a generic quark 
doublet {U,D). For the U quark, we have 



Smu 
mu 



Xu 

sin^/? \2 

3 

Xu ' 



^A + ^lnAr + 4 



m,r 



+ 



sin /? cos^ /3 
^^ ^^A + iln^ 



^ cos^ /3 



-A 
) 

3 

+ 4 



1, /"' 
- in — - 

2 m' 



+ 



u 



if mu > mn, 
if m^u = TTiD, 
if mu <C mD- 



(34) 
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The corresponding expression for the D quark is obtained from Eq. (1M|) by interchanging 

iriu <-i> mo and sin (3 *-^ cos f3. Furthermore, we have 



fNr 



-Xu 



if mu > mo, 
if m^u = moi 



(35) 



Nr 



xd if niu -C mo- 



The counterparts of Eqs. fl34l) and fl35l) for a generic lepton doublet (A^, E) are obtained 
by substituting m^u — >■ ^at, '^d — ^ ^e, and A^c — >■ 1- 



Our final result for the three-generation case reads 



^r = 2^/^Gf ^iV.x. cot /5 (^ + ^ - 1^^ - i^^ + ^ cot^ 



6 



TT 



„ „ + — cot^ /5 + - cot^ 3 

9 9 sin^/? 9 sin^/? 9 ^ 9 ^ 



9 ' 



2l/4^1/2«em^^ 

vr 



Xj 



14 2 

yCOt=^ /? + -(? -AT,) cot/? 



(36) 

(37) 



where a is the angle that rotates the weak eigenstates of the neutral CP-even Higgs bosons 
into their mass eigenstates. On the right-hand side of Eq. (I36l) . we exhibit separately the 
finite contributions from the x^i 0^? h^ i H^i ^^^ i and H^ bosons, and the universal 
counterterm, after top-quark mass renormalisation, in the order in which they appear 
there. The UV-divergent parts vanish after exploiting simple trigonometric identities. 
The same is true for the unrenormalised contributions from neutral particles as well as 
for the naive contributions from the asymptotic expansion of the diagrams containing 
charged particles. However, the total result is nonzero, as opposed to the QCD case. We 
observe that the a dependence carried by the contributions from the neutral CP-even 
Higgs bosons cancels in their sum. This reflects the fact that, by neglecting their masses, 
we effectively treat the h^ and H*^ bosons as mass degenerate, so that we may rotate the 
angle a away. 

Now we examine the influence of a sequential generation consisting of heavy fermions 
F as specified above. Besides the appearance of new leading correction terms quadratic in 
their masses, of generic order 0{xf)i also the 0{xt) correction is then modified. This may 
be understood by observing that the LO result then receives three more mass-independent 
contributions in addition to the one from the top quark, from the charged fermions f/, D, E, 
which feed into the 0{xt) correction through the universal counterterm 6u- This may be 
accommodated in Eq. (!36l) through the substitution 



9 ' 



2 2 

9 9 



— - tan^ B tan^ B. 

18 '^ 2 ^ 



(3^ 
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The 0{xf) contribution due to the (f/, D) doublet is found to be 



{U,D) 



21/46-^ (a,^/vr)iV, 



,4 3 ^ 13 - 4Ar, ^ 7 + A^e 

^u I o ^o* P "^ f; — ^o^ -^ 73 — ^^^P 

6 y io ^ 

4 o 17 8 1 ^ 

X{/ I - cot /? H cot /5 H — tan {3 -\ — tan {3 



9 



xd -{I -N^) cot (3 + 



9 



9 

tanp H — tan p 
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1 
3 



iimu > m^i, 
if mjj = m£,, 

if m(/ -C m£). 

(39) 

In each case, we find that the contributions from the various proper two-loop diagrams 
cancel, so that we are only left with the counterterm contributions. Again, the total result 
is non-zero, in contrast to the QCD case. 

The counterpart of Eq. fl39|) for the (A^, E) doublet is simply obtained by appropriately 
adjusting the quantum numbers Nj and Qf and reads: 



AXF 



2y^G';\a,Jn) 



= < 



Xn 
Xn 
Xe 



l-^ArJcot/3+i 



1 ] tan/3 

9 ' ^ 



(cot (3 + Atan(3 + 3 tan^ /3) 



4 1 / A^^ \ 
— NcCotp + - I 5 -] tan/3 + 3tan^/3 



if m-AT > m^, 
if rriN = mE, 
if mjv <^ ruE- 
(40) 



3.3 0{xt) and 0{xf) corrections to r(A° -^ gg) 

Now we turn to the 0{xt) and 0{xe) corrections to r{A^ -^ gg). They may be easily 
extracted from the analogous calculation for r(y4° -^ 77) discussed in Section 13. 2[ by 
retaining only those diagrams where both photons couple to quark lines and substituting 



"em^; g^; g'-' '( 



a. 



Note that this does not affect the factors of N^ originating from the 



renormalisation procedure. 

In the three-generation case, we thus obtain 



AT 



2'/^GF—Xt cot /3 f 5 + 3 - 3^^ - 3^^ + 5 cot^ /S + 3 cot^ /? - ^ ) (41) 
Tc \ sin (3 sin (3 2 



2^/'GE-Xt 



5 cot^ /5 + 5 - 



N. 



cot/? 



(42) 



where the five terms on the right-hand side of Eq. ( HTi) again represent the finite contribu- 
tions from the x°) 0^5 ^"5 H^, ^°5 and H^ bosons, and the universal counterterm, after 
top-quark mass renormalisation. The final result is again independent of a. 

In the presence of a sequential heavy-fermion generation, the universal counterterm in 
Eq. (14T|) is modified according to 



2 ~* 2 2 2 ^^ ^' 



(43) 
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The 0{xf) contribution due to the (f/, D) doublet is found to be 



AXF 



Xu 



3 cof' /3 + (4 - AT,) cot /3 + 1 



Nr 



tan/5 



2V4G]/^(a,/7r) 



Xu (3 cot^ /3 + 5 cot /3 + 5 tan /3 + 3 tan^ (3) 



Xd 



iVc 



;i - A/'c) cot /5 + ( 4 - — ) tan/3 + 3 tan^ /3 



if mu > rriD, 
if mu = mj:,, 
a mu <^ mi). 



(44) 

The {N, E) doublet can generate a 0{xf) contribution only through the universal 
counterterm, so that 



AXF 



2V4Gf(a,/7r) 



xn ( — cot 13 — - tan (3 j if m,j<i ^ tue, 

if m^v = ■?«£;, 
xe ( — cot (3 — - tan /5 J if m^v -C mg. 



(45) 



Note that the 0(q;s) correction to r(y4'^ -^ gg) cannot be recovered from the one to 
r(y4° — > 77) because it receives additional contributions from diagrams involving gluon 
self-couplings. 



4 0{xf) correction to T{H — ^ 77) 

Applying similar techniques as in Section [221 we now also derive the 0{xf) correction to 
r{H — >■ 77) in the SM endowed with a sequential generation of heavy fermions. Due to 
electromagnetic gauge invariance, the transition-matrix element of if ^ 77 possesses the 
structure 

r = [(gi • q2)g'"' - qU^K(qi)e:(q2)A. (46) 

To obtain a strong check on our analysis, we actually verify electromagnetic gauge in- 
variance by separately projecting out the coefficients of the Lorentz tensors (gi ■ q2)g^'^ 
and q^q2 in Eq. (HBl) . Furthermore, we work in general R^ gauge, so as to verify that the 
gauge parameter ^ cancels in the final result. From Eq. fH^ . we obtain 



TiH 



77) 



Ml, ,,n 



(47) 



where Mh is the mass of the SM Higgs boson. 

The form factor A is evaluated in perturbation theory as 

A = A]^ + Y, {-^f + Af+Ay + ■■■)+■ 
f 



(4J 



where A}^ denotes the one-loop contribution due to the W^ boson and the other contri- 
butions carry similar meanings as in Eq. (I23ll . 
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A comprehensive review of the present theoretical knowledge of r{H -^ 77) may be 
found in Ref. [lOj. The LO result was first obtained in Ref. [H]. The 0(a,) [lilli2] and 
0{a1) [43] QCD corrections are also available. As for the two-loop electroweak correction, 
the contributions induced by light [44J and heavy fermions [20|45] as well as the residual 
ones [16] were recently evaluated. The 0{xf) correction due to a sequential generation of 



heavy fermions was studied in Ref. pO] for general values of their masses. In the following, 
we revisit this analysis for the mass hierarchies rriu ^ mo, TTtjj = m£,, and mu <^ rriD, 
and similarly for the leptons A^ and E using asymptotic expansion. We refrain from 
considering r{H -^ gg). 
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Figure 6: Generic two-loop diagrams contributing to the 0{xf) electroweak correction to T{H — > 
77). S = x*^, (t^^ , W^, H and F = U,D,N,E denote generic bosons and fermions, respectively. 
The couplings of the neutral particles to the photon vanish. 



The diagrams contributing to r{H -^ 77) at 0{xf) are shown generically in Fig. [6l 
In contrast to the case of A° -^ 77, now also virtual W^ bosons participate in 0{xf) 

mm. 



Unlike in Refs. [TTPS] , we choose not to perform tadpole renormalisation here. In 
turn, we need to include the diagrams that are generated by attaching a Higgs tadpole 
in all possible ways to any one-loop seed diagram. Some examples are depicted in Fig. [3 
These tadpole diagrams yield terms proportional to nip, which cancel against the tad- 
pole contributions to the counterterms, which are all proportional to m^, and the nip 



16 




H 



1 





H 



F)S 



H 



S 



S 



7 



7 




Figure 7: Examples of two-loop tadpole diagrams contributing to the 0{xf) electroweak cor- 
rection to T{H — > 77). S = x*^, 'Z'^) W^^j -ff and F = U,D,N,E denote generic bosons and 
fermions, respectively. The couplings of the neutral particles to the photon vanish. 



terms from the Higgs mass renormalisation and the asymptotic expansion of the two-loop 
diagrams involving virtual 0^ bosons. Thus, the final result is devoid of m"^ terms. 

It is interesting to note that the contribution from the attachment of a tadpole to 
a fermion line is cancelled by the nip term from the renormalisation of the mass of 
that line in the one-loop seed diagram. Furthermore, the tadpole contributions to the 
renormalisations of the factors 1/Mw and mp in the HFF vertex cancel each other. 
These are the only tadpole contributions that could generate m?p terms in the final result, 
through the expansion in Tp- 

Prior to evaluating the proper diagrams of Figs. [6] and [71 we discuss the renormalisation 
in some detail. As before, we need to renormalise the Higgs-boson wave function and the 
masses of the W^ boson and the heavy fermions. In addition, we now also need to 
renormalise the Higgs-boson mass. The parameter M\y appears in the W^, (p^, and 
u^ propagators and in the HW^W^, Hcp^W^, (p'^W'^'y, Hcp'^W'^'j, and HFF vertices, 
where u^ are the charged Faddeev-Popov ghosts. The only vertex involving Mh is Hcp'^cfy^ , 
which induces two-loop contributions via 6Mfj- Finally, rrip occurs in the F-fermion 
propagator and in the HFF vertex. The corresponding counterterms are defined through 



ITLp 



mp + 5m p + dm^p^, 

Ml + 5 Ml + ^M^**"^, 
ZhH = ( 1 + Uzh 



H, 



(49) 
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where tadpole contributions are marked by the superscript "tad." Note that 5Zh is 
obtained from the derivative of the Higgs-boson self-energy and thus has no tadpole 
contribution. 

In the case of the (f/, D) doublet, the counterterms in Eq. (H9l) read: 

,2 \ 

iimu > rriD-, 



Sm 



u 



^'''^^ + l'°£!;^' 



mu 



2xu 
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Xd 
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-1- 


2 mi, 


5 
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if mu = moi 



5m^^'^ 
mu 



2 / 2 

xc/A^ctS- I 4A + 4 In -^ + 4 ) if m[/>m^, 



M: 



¥ 



xuN, 



m 



"^"^^8A + 81n^ 
mf 



M| 



V 



+ 8 if mu = rriD, 



m 



D 



^H 
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4 A + 4 In -^ + 4 iimu <^ mo 
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iimu ^ f^D, 
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-24A-241n-^-8 ) +4A + 41nAr-- 
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'^ / mi, 3 



m 



D 



D 



if mu > m^, 
if mu = m,D, 
a mu <C mo, 



SZ 



H 



xuNc{-2A-2ln-^ + -] if mc; > m^^, 
mp 3. 

a^t/A^'c ( -4A - 4 In -^ + ^ ] ii mu = mo, 
mfj 3 

u^ 4' 

xdN^ [ -2A - 2 In -V + o ) if ^(7 < m^ 
mj-, 3 



(50) 



and similarly for dm^D/m^r, and 5m^^'^ /mr,. The renormalisations of the Higgs-boson wave 
function and the IV-boson mass in the HFF Yukawa coupling combine to a universal 
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correction 



(^u = 2 ( ^^H 



5M^ 



(51) 



which should be compared with Eq. fl5U]) for the case of the A^ boson. We have [SO] 

(7 



5u 



xuNc a niu 3> mn, 



xuNc if mu = mo, 



6 

4 

3 

7 





(52) 



The counterparts of Eqs. (ISUj) and (132]) for the (A^, E) doublet are obtained by substituting 
rriu — > ttitv, '"^d -^ tt^Ei and iVc ^ 1. Those for the three-generation SM may be found in 
Eqs. (68)-(71), (74)-(76), and (80) of Ref. [i8|. 

We now list our final results for J^/ due to the {U, D) doublet. For the sake of 
comparison with Ref. [20] , we exhibit the dependence on the electric charge of the heavier 
quark, exploiting the relation Qjj = Qo + 1. We find 



AXF 

{U,D) 



2^/^Gf{a,j7i)N, 




2^ 7 

-:^-QQu + ^Ql - ^N, (1 - 2Qu + 3Ql) 

-^ - 8Qu + 8Ql - ^N, (1 - 2Qu + SQfj) 

2^ 7 

-gg + 6(3d + iQl - -jN, (2 + iQo + 3Ql) 



12 9 



9 

Nc 

4 9 



if mu = rno-i 
iimu <^ rriD- 



if mjj = rrij:,, 
a mu <C mu, 

(53) 



Appropriately adjusting the quantum numbers in the various contributions, we find 
Ay due to the (iV, E) doublet to be 



AXF 

•^{N,E) 



2y^Gf{a,Jn) 




if niN > niE, 
if rriN- = niE, 
if rriN -C niE. 



(54) 



In the remainder of this section, we compare our results with those obtained in 
Ref. [20]. That reference is more general than ours in the sense that no hierarchies 
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among the heavy-ferinion masses are assumed. However, we can compare Eq. (!53l) with 
Eqs. (56) and (57) of Ref. [20], where the hierarchies rriu 3> m_D and rriu = m£,, respec- 
tively, are considered. We can reproduce these equations if we include an overall minus 
sign in our expression for 6mu/mu in Eq. (l50l) . In other words, there should be an overall 
minus sign on the right-hand side of Eq. (9) in Ref. [20], which was already noticed in 
Ref. [l5] in connection with the SM case. In Ref. ^5j, it was also observed that the limit 
■m£, — > of the fourth-generation result for M^y -C mo <^ rriu differs from the calcu- 
lation with mo = from the beginning, which is appropriate for the third generation, 
where rrib <C Mw <^ m,t. These two observations lead to a modification of Eq. (60) in 
Ref. [20], where the leading correction due to the five heavy fermions t, U, D, N, and E 



is presented. Adopting the notation of Ref. [20], we have 
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(55) 



where a 



"U ""D ""D, 

aem and G^ = Gp. By the same token, Eq. (59) of Ref. [20] becomes [I5|l6] : 



A: 
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aGl^^ 47 



7r 



23/4 9 



1 + 
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M 



87T^V2 



367 



94 



-m. 



(56) 



In Ref. [20], the 0{xf) corrections to the Hgg coupling in the SM with and without 
a sequential generation of heavy fermions were inferred from those if 77 diagrams were 
the photons are directly coupled to loop quarks. In the case of the 0{xt) corrections, the 
effects due to the flipped sign in the top-quark mass counterterm and the interchange of 
mass limits in the proper diagrams incidentally compensate each other, so that Eq. (61) 
in Ref. [20] agrees with Refs. [l5ll5T] . In the case of the 0{xf) corrections, where no mass 



limits are interchanged, the fermion mass counterterms cancel within each isodoublet, so 
that Eq. (62) of Ref. [2U] goes unchanged. 



5 Discussion 

We now explore the phenomenological implications of our results for r{A^ -^ 77) and 
r{A^ -^ gg). For definiteness, we concentrate on the more likely case of three generations. 
As explained in Section [1], we consider a scenario with low to intermediate values of the 
Higgs-boson masses and tan/3 and large values of the supersymmetric-particle masses, 
so that the dominant electroweak two-loop corrections are of relative order 0{xt). We 
adopt the following values for our input parameters [31j: Gp = 1.166 37 x 10^^ GeV~^, 
as{Mz) = 0.1176, and m,t = 170.9 GeV. As for the unknown 2HDM input parameters, 
we assume that M40 < 160 GeV and 2 < tan/3 < 10. For larger values of tan/?, our 
approximation of neglecting the bottom-quark contributions is likely to break down. 
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Figure 8: 0{xt) and 0{as) corrections to r(^'' — > 77) (a) for M^o = 100 GeV as functions of 
tan P and (b) for tan /? = 2 as functions of M^o . 
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77). Neglecting the bottom-quark contribution, its 0(xt) 



4 + 



tair (3 



Oin) 



(57) 



where, in the second equality, Al''~' is approximated by the leading term in the second 
line of Eq. (12^ for f = t. We observe that this correction is negative, has its maximum 
size for small values of tan/3, and is independent of M^o, apart from the Mj^o dependence 
carried by ^J"°. Its evaluation according to the first line of Eq. fl57p is compared with 
the 0{as) correction in Fig. [HI We observe from Fig. M^a) that the 0{xt) correction 
amounts to —1.7% at tan/5 = 2 and rapidly reaches its asymptotic value of —1.2% as tan/3 
increases, whereas the 0{as) correction, evaluated from Eqs. (^^ and fl27j) . is positive 
and independent of tan/3, as long as the bottom-quark contribution is neglected. The 
Mao dependence of the 0{xt) correction shown in Fig. Mh), which is induced by A]"^ 
as mentioned above, is rather feeble, so that we may expect the unknown 0(r") {n = 
1, 2,3,.. .) terms in Eq. (1571) to be of moderate size, too. The smallness and approximately 
quadratic M^o dependence of the 0{as) correction is due to the absence of the leading 
0(r°) term in A"", given by Eq. (127|) . We conclude that the 0{xt) reduction more than 
compensates the 0{as) enhancement for M40 ^ 120 GeV. 

We now turn to T{A^ — >• gg). For Tj <^ 1, its 0(xj) correction reads 
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As in the case of ^4° — ^ 77, this correction is negative, has its maximum size for small 
values of tan/?, and is independent of M^o, apart from the M^o dependence carried by 
A}''"'. In Fig. [9l its evaluation according to the first line of Eq. (!58l) is compared with the 
full 0{as) correction [10] due to virtual top quarks, which also involves three-parton final 
states. In contrast to the case of A^ — > 77, the 0{as) correction to F(A° -^ gg) does 
have a 0{t^) term, which is about 68%. The 0{t^) term is also known at 0{al), where 
it is still as large as 23% fTB]. The 0{xt) correction to r(A'' — >■ gg) ranges from —2.8% 
at tan/? = 2 to the asymptotic value —2.1% and partly screens the sizeable 0{as) and 
0{a1) corrections. 

Let us now briefly comment on the 0{xf) corrections to F(y4° — > 77) and F(y4° —>■ gg) 
due to a sequential generation of heavy fermions, given in Eqs. f p9|l . (HOj) . (jH]), and fj45l) . 
These can be sizeable for large values of mp just because of the prefactor xp- For large 
values of tan/?, further enhancement comes from the terms of maximum power in tan/?, 
which are cubic for mu ^ rriD, fnu <^ mo-, ttin ^ itie, and rriN <^ tue in the case of 
A^ — >• 77 and for mu ^ rriD and mu <^ m,£) in the case of A^ -^ gg, and (at most) linear 
for the other mass hierarchies. For tan/? > 2, the 0{xf) corrections reduce the LO results 
for r(A° — s> 77) and T{A^ -^ gg), except for the corrections due to a {U, D) doublet with 
m,u ^ rrij^ and, in the case of V{A^ -^ gg), also those due to the (iV, E) doublet. Because 
of the constraint from electroweak precision tests [31] on the rho parameter [20|52] , the 
case of approximate mass degeneracy within the {U, D) and (A^, E) doublets is favoured, 
so that a tan^ /?-enhanced screening is likely to be encountered. 



6 Conclusions 

In conclusion, we analytically calculated the dominant electroweak two-loop corrections, 
of order 0{xt), to r(A° —^ 77), T{A^ -^ gg), cr(77 -» A^), and a{gg -^ A°) within the 
2IIDM with low- to intermediate- mass Higgs bosons for small to moderate value of tan jS 
using asymptotic expansion in M^o/(2mi)^. We also studied how these corrections are 



22 



modified by the presence of a sequential generation of lieavy fermions, with generic mass 
mp, and provided the 0{xf) corrections arising then in addition. We also revisited the 
0{xt) and 0{xf) corrections to T{H — *• 77) and (7(77 — > H) in the four-generation SM 
and clarified an inconsistency in Ref. [20] . 

We recovered the notion that the naive treatment of the 75 matrix being anticom- 
muting in d space-time dimensions leads to ambiguous results, which depend on the 
way of executing the Dirac traces. To consistently overcome the non-trivial 75 problem 
of dimensional regularisation, we adopted the HVBM scheme [2T] and included a finite 
renormalisation constant, Z^, for the pseudoscalar current to effectively restore the anti- 
commutativity of the 75 matrix [22|23f24] . The 0{xt) and 0{xf) terms of Z^ were found 
to vanish. We worked in the electroweak on-shell renormalisation scheme [321133] endowed 
with the DCPR definition of tan/3 [351136] . 

On the phenomenological side, the 0{xt) correction to r{A^ — >■ 77) and (7(77 — *> A°) 
is of relative importance, since it more than compensates the 0{as) enhancement for 
M^o ^120 GeV. It leads to a reduction of the LO results, which ranges between —1.7% 
and —1.2% for 2 < tan/3 < 10 and is independent of M^o. Such an effect might be 
measurable for (1(77 — > A^) at the ILC operated in the 77 mode on the A^-boson resonance 
and possibly also for r{A^ —>■ 77) at the ILC in the regular e^e~ mode [7]. 

As for r{A^ -^ gg) and (j{gg — > A"), the 0{xt) correction screens the sizeable QCD 
enhancement, by between —2.8% and —2.1% for 2 < tan/5 < 10, and is independent 
of Mj^o. Such a reduction of cr{gg — » A^) should matter at the high luminosities to be 
achieved at the LHC E.g., given an annual luminosity of 100 fb~^ per LHC experiment, 
a. pp ^ A^ + X cross section of about 35 pb (for M^o = 100 GeV and tan/3 = 2) [3] 
amounts to 7 x 10^ A" bosons per year, 2.8% of which still corresponds to a substantial 
subsample of 200.000 A^ bosons per year. Furthermore, the size of this correction is in 
the ballpark of the theoretical uncertainty due to the parton distribution functions [ID] 
and the scale dependence of the NNLO QCD prediction [T7] . 
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